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' We study two-component solitons and their symmetry-breaking bifurcations (SBBs) in linearly 

coupled photonic systems with a spatially inhomogeneous strength of the coupling. One system mod- 
els an inverted virtual photonic crystal, built by periodically doping the host medium with atoms im- 
plementing the electromagnetically induced transparency (EIT). In this system, two soliton- forming 
probe beams with different carrier frequencies are mutually coupled by the EIT-induced effective 

■ linear interconversion. The system is described by coupled nonlinear Schrodinger (NLS) equations 
for the probes, with the linear-coupling constant periodically modulated in space according to the 
density distribution of the active atoms. The type of the SBB changes from sub- to supercritical 

■ with the increase of the total power of the probe beams, which does not occur in systems with con- 
stant linear-coupling constants. Qualitatively similar results for the SBB of two-component solitons 
are obtained, in an exact analytical form, in the model of a fused dual-core waveguide, with the 
linear coupling concentrated at a point. 
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The effect of the spontaneous symmetry breaking of solitons in two-component linearly-coupled system was studied 
in detail in nonlinear optics [l[ and Bose-Einstein condensates (BECs) 0, HJ ■ The symmetry-breaking bifurcation 
(SBB) occurs as a transition of a symmetric localized (solitonic) ground states into an asymmetric one when the 
linear-coupling constant drops below a critical value. In optics, the linear coupling originates from the overlappingof 
evanescent fields between adjacent waveguides, such as in dual-core fibers pj, HTHI and arrays of such fibers [E 0], 
Q\ • or from the linear mixing of orthogonal polarizations induced by the twist or elliptic deformation in bimodal fibers 
@, 0. In binary BECs, a similar effect originates from the interconversion between hyperfine atomic states induced 
by a resonant electromagnetic wave, as demonstrated theoretically in a variety of settings floj |. There are two kinds of 
| the SBBs, sub- and supercritical. In the case of the subcritical symmetry breaking (which is tantamount to the phase 
f^i . transition of the first kind), the system features branches of asymmetric states which emerge as unstable ones, going 
04 ' at first backward from the bifurcation point and undergoing the stabilization after turning forward. In the case of 
, the supercritical SBB (it is tantamount to the phase transition of the second kind), asymmetric branches emerge as 
stable ones and immediately go in the forward direction. An interesting problem is a possibility to control the type of 
the symmetry breaking, and thus to switch between the respective phase transitions of the two kinds. For example, 
the addition of a periodic potential (optical lattice) acting in the unconfincd direction changes the character of the 
SBB from sub- to supercritical [H), and a similar effect is induced by rendering interactions nonlocal 11|. 



One of versatile techniques for the control of transmission properties of optical media is the electromagnetically 
induced transparency (EIT) [l2| . The EIT gives rise to a variety of nonlinear features, which can be used for the 
making of single- and multi-component solitons, and thus, subsequently, for the implementation of the symmetry 
breaking in solitons. These features include the self-enhanced [l3| and giant [l4| Kerr effects, as well as the enhanced 
frequency conversion [l5l - [T8j . 

The first aim of the present work is to study the SBB for solitons in two-component photonic systems featuring 
spatial modulations of the strength of the linear coupling induced by the EIT-mcdiatcd frequency conversion. One 
such system represents a virtual photonic crystal (PhC) formed by a periodic modulation of the concentration of 
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(a) (b) 

FIG. 1: (Color online) (a) The energy-level structure (the double-A scheme) of the active dopants (Pr 3+ ). (b) The periodic 
distribution of the dopants (filling blue regions) against the background (gray regions). The dopant density is represented by 
function R(x) 



active atoms doping a passive host medium. This technique has been recently implemented in the fabrication of 
imaginary-part PhCs by implanting atoms of RhB (Rhodamine B, a dye manifesting saturable absorption) into the 
SU-8 polymer, which is a commonly used transparent negative photoresist (T^ |. The present model refers to another 
pair of the active atoms and passive background, namely, Pr 3+ ions and the YSO crystal, respectively. 

The energy scheme of Pr 3+ is displayed in Fig. [TJa), which is a four- level energy system. Fields Hi and f^, which 
have different carrier frequencies, represent two soliton-forming beams. They induce resonant transitions from the 
same ground state, i.e., |1), to different excited states, |3) and |4), with detuning A x and A 2 , respectively. The 
resonant transitions from a common metastable state, |2), to the same pair of the excited states (|3) and |4)) are 
powered by pump (coupling) fields Clci and fic2, with detunings adjusted so that Aci = Ai and Ac2 = A 2 , which 
nullifies the two-photon detuning. The four probe and pump fields involved into the scheme build a double-A scheme 
fl8l |. In this setting, the transitions driven by coupling fields flci and Oc2_give rise to the effective linear mixing 
between fields ip an d 0, in the form of the induced frequency conversion (l7l . Il8j . The latter feature places the system 
into the class of the linearly-coupled two-component ones. 

The periodic modulation of the dopant concentration, i.e., the structure function R(x\ which represents the density 
of the implanted ions in Fig. [ljb), turns the uniform medium into a virtual PhC [20|-(23l|. so called because its 
properties are controlled by the pump fields, rather than by a permanent material structure (other well-known 
examples of virtual structures are lattices optically induced in photorefractive crystals [24j|). The density- modulation 
pattern induces an effective periodic linear potential, along with a spatially periodic modulation of all local optical 
parameters, including the strength of the linear interconversion and coefficients of the self-phase modulation (SPM). 
The virtual PhC of this type may be naturally called an inverted one, in comparison with ordinary (material) PhCs 
built in self-focusing media: While the ordinary PhC structures feature effective linear and nonlinear potentials with 
coinciding local minima, in the present setting minima of the local potential coincide with maxima of its nonlinear 
counterpart, and vice versa, i.e., the nonlinear potential is inverted with respect to the linear one (2pj-[22l|. Similar 
models with competing effective linear and nonlinear potentials were introduced, in a different context, in Ref. [25| . In 
fact, this setting may be considered as a variety of the general concept of nonlinear lattices and mixed linear-nonlinear 
ones [26[. 

The studies of the single-component model of the inverted virtual PhC with the Kerr and saturable nonlincaritics, 
reported in Refs. [2(| and [2l[, respectively, have revealed that the stable location of the solitons in such systems is 
controlled by the total power, with the low- and high-power solitons tending to be pinned by minima of the linear 
and nonlinear potentials, respectively. Transitions between the different positions of the solitons were identified as 
supercritical SBBs. 

The major objective of this work is to study two-component solitons, supported by the EIT scheme, in the inverted 
crystals, the main point being the SBBs of the solitons. In Sec. II, we derive the model by means of the semi-classical 
consideration of the interaction of the electromagnetic fields with atoms. Then, we focus on effects of the periodically 
modulated linear coupling between the two components, which is the main novel element of the model. In Sec. 
Ill, we consider the competition of the linear coupling and SPM nonlinearity. We concentrate on the SBB specific 
to the two-component system, and do not dwell on the above-mentioned breaking of the spatial symmetry of the 
solitons, which can be adequately studied in the single-component system [2lL 1221] . A systematic numerical analysis 
demonstrates a transition of the subcritical SBB into a supercritical one with the increase of the soliton's total power, 
and a concomitant shrinkage of the soliton. In Sec. IV, analytical results are presented for an allied model of a 



3 



dual-core fused [27[ spatial coupler, with the tightly concentrated linear coupling represented by the delta-function 
(in fact, the same system, but with a periodic array of coupled sites, and a periodic modulation of the refractive index 
and Kerr coefficient, may also represent the model of the virtual PhC considered in Sec. II). In that model, the SBB 
of two-component solitons can be studied in an exact analytical form. The paper is concluded by Sec. V. 



II. THE MODEL OF THE INVERTED NONLINEAR PHOTONIC CRYSTAL 



is 



The atomic Hamiltonian corresponding to the configuration shown in Fig. 1(a) 

H = &[<S|2)<2| + Ai|3><3| + Aa|4><4|] 

-a[fii|3)<l| + n 2 |4)<l| + n C i|3)<2| + fl C2 |4)(2| + H.c], (1) 

where fields Q and detunings A 12 are defined as per the figure, H.c. stands for the Hcrmitian-conjugate contribution, 
and the two-photon detuning, S, is set to be zero. Under physically realistic conditions, spontaneous decay of the 
states included into the scheme may be neglected. Further, if the dopant atoms are initially kept in the ground state, 
steady-state solutions for density-matrix elements, p^i and P41, which are activated by the two soliton-forming probes, 
can be written as follows: 

a 2 |^| 2 |o 2 |2 \n C2 f n cl ^ 2 n 



Ai|Q 2 | 2 |«i| 2 |n C i| 2 gcgQ^ 

a L(~t i O C C 



where fl^ = A 2 |f2c?i| 2 + Ai|f2c2| 2 - The first (SPM) terms on the right-hand sides of Eqs. originate from the 
self-enhanced Kerr effect [llj], the second and the third (XPM) terms represent the giant Kerr effect [l4j . and the 
fourth terms represent the linear coupling between the soliton-forming beams, which originate from the EIT-induced 
frequency conversion fl7| . 

The polarization experienced by the two probes in the medium are [28j 

& 1 (x)=2N{x)p 3lPsi , 

3* 2 {x) = 2iV(x)p 4 iP4i. (3) 

Here,p3i and p4i (which are assumed real) are the matrix elements of the dipole transitions |1) — > |3) and |1) — > |4). 
A reasonable simplification of the model is attained by assuming that p3i w pn = p. The paraxial propagation 
equations for the slowly varying envelopes of the probe fields are 

Substituting Eq. ^ and Eq. ([3]) into Eqs. ((3]), one arrives at coupled nonlinear Schrodinger (NLS) equations, 

= -ld xx ip + R(x) (Vnp + aMH' + + Ci4>) , 

ldz(j) = -ld xx <j> + R(x) (U 2 + a 2 |0| 2 + n\^\ 2 4> + C 2 V) • (5) 

where = ^4/7, - fi 2 / 7 , R(x) = p 2 N (x) / e h~/ 2 , V 1 = -\n C2 \ 2 /Q 2 c , V 2 = -|fici| 2 / n c. tr x = A 2 /A 1 Q 2 C , 
(j 2 = Ai/A 2 f2 2 7 , and k = l/tig. If we let fici and ^C2 be real, then coefficients C\ — C 2 = C = Qci^C2/^c 
account for the EIT-induced linear mixing of the probe fields. The effective linear-coupling coefficient CR(x) in Eq. 
([5]), which is periodically modulated due to the distribution of the dopant density, makes the system different from 
various previously studied models of linearly-coupled systems @- 



III. NUMERICAL RESULTS FOR THE MODEL OF THE PHOTONIC CRYSTAL 



To focus on effects of the periodically modulated linear coupling competing with the SPM terms, we drop the XPM 
interaction in Eq. (5), and simplify Eqs. ([5]) to the following form: 

id x i> = -\d xx i> + V(x) (1 - |^| 2 ) i, - C(x)<f>, 

id z <t> = ~d xx cf> + V(x) (1 - |0| 2 ) - C{x)4>. (6) 
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FIG. 2: (Color online) Stationary profiles of the two components of stable symmetric (a) and asymmetric (b) solitons, both 
found for P = 4, Co = 0.078. These solitons belong to the bifurcation diagram displayed in Fig. HJa). In this figure and below, 
the black solid lines display the shape of the underlying modulation function, R(x) [see Fig. QIb)]. 
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FIG. 3: (Color online) Examples of stable symmetric (a) and asymmetric (b) solitons found for P = 10 and Co = 2 (a) or 
C = 1.5 (b). These solitons belong to the bifurcation diagram displayed in Fig. 2£d). 



These equations generalize those derived in Refs. [2(| Hlj for the inverted PhC with the 7r-shift between the periodic 
linear and nonlinear potentials, under the assumption that the depletion of the coupling fields, flci and f2c2 [see 
Fig. HJa)], may be neglected. In this work, we consider the density-distribution profile in Fig. [TJb) corresponding to 
R(x) = cos 2 x, i.e., V(x) = Vq cos 2 x, C(x) — Co cos 2 x, the notation being fixed by scaling the modulation period to 
be 7T. Generic results are displayed below for amplitude Vq = 0.5 of the nonlinearity modulation, while strength Co 
of the linear coupling is varied. As for the sign of Co, it may be fixed to be positive, as Co < can be transformed 
into Co > by the change of <j> — > —(f), while the sign of ip is not altered. 

Stationary soliton solutions to Eqs. (|6|) were found in a numerical form by means of the imaginary-time propagation 
method [29[, and their stability was subsequently tested by direct simulations of the perturbed evolution in real time. 
The solitons are characterized by the total power, 

/+oo 
[|V(x)| 2 + \4>{x)\ 2 } dx. (7) 
-oo 

As expected, both symmetric and asymmetric soliton modes, in terms of the coupled components, were found, see 
Figs. [5]and[3] These examples display stable solitons, which are broad in comparison with the underlying modulation 
pattern at smaller values of P [Fig. [2], and narrower modes, with the width comparable to the modulation period, at 
larger P [Fig. [3]. 

The numerical results are summarized in the form of bifurcation diagrams displayed in Fig. U at different fixed 
values of the total power, P. The bifurcations are driven by the decrease of the coupling coupling, Co- Unstable 
branches of symmetric solitons, which should continue the stable ones in all the panels, and narrow intermediate 
branches of unstable asymmetric solitons in panel Eta), are missing as the imaginary-time integration method does 
not converge to unstable solutions. The diagrams show the transition from the subcritical SBB to the supercritical 
bifurcation with the increase of P. 
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FIG. 4: (Color online) The bifurcation diagrams for the two-component solitons with different fixed values of the total power, 
corresponding to the variation of strength Co of the linear coupling. Panels (a,b) and (c,d) correspond to the bifurcations of 
the sub- and super-critical types, respectively. 



For broad solitons one may replace R(x) = cos 2 x in Eqs. ([6]) by its mean value, (cos 2 x) = 1/2, which reduces the 
system to the one studied in earlier works [l|, where the SBB is subcritical, in accordance with Fig. HJa). The new 
situation is actually found here for narrow solitons, for which the SBB turns out to be supercritical. 

It is instructive to compare these results with those found in a system where the linear and nonlinear potentials 
are the same as in Eqs. (JH]), but the linear coupling is made constant, by replacing Cocos 2 (a;) — > Co/2: 

id z i> = -\d xx i> + v(x) (i - m 2 ) V - \c Q <t>, 

id z ct> = -\d xx <\> + v(x) (i - H 2 ) <i> - ^Coi>. (8) 

Figure [S] demonstrates that the modified system docs not give rise to the transition of the subcritical SBB into the 
supercritical type, even in the case when P is large and the solitons are narrow. Thus, the periodic modulation of the 
linear coupling is essential for this transition. 



IV. ANALYTICAL RESULTS FOR THE FUSED COUPLER 



Numerical finding presented in the previous section suggest that the spatial modulation of the linear coupling is a 
crucial factor which determines the kind of the phase transition (SBB) in the system, and the possible change of the 
kind. In this section, we aim to illustrate the genericity of this rule by means of exact analytical results, obtained in 
an allied system which describes a fused dual-core nonlinear planar waveguide [27| , with the linear coupling tightly 
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FIG. 5: (Color online) The bifurcation diagrams for the two-component solitons produced by modified system ((8]). 



concentrated at x = 0: 



1 2 

id z ip = --d xx %l)-\il)\ ^-28{x)4>, 

1 2 

id z (j> = --d xx <j>-\4>\ 4>-28{x)^, 



(9) 



where 8(x) is the delta-function. This system is a limit case of Eqs. ©, with C(x) = 28{x) and V(x) = 1 (in that 
case, the constant linear potential is trivial in the present case and may be dropped). Coefficients in the system can 
be fixed as in Eqs. ^ by means of an obvious rescaling. It is relevant to mention that symmetric and asymmetric 
solitons in a discrete version of this system were recently studied in Ref . , but the discrete system does not admit 
exact solutions. In fact, the underlying system of Eqs. © may also be realized in terms of the dual-core spatial 
coupler, with periodically modulated strength of the coupling, refractive index, and Kerr coefficient. 

Stationary solutions to Eqs. © are sought for as {ip, (j)} = ex P (ikz) {u(x), v(x)} , with real functions u(x) and v(x) 
satisfying equations 



ku 
kv 



-28(x)v, 
28{x)u. 



(10) 



Obviously, solutions to Eqs. (| 10[) are subject to the following boundary conditions, produced by the integration in an 
infinitesimal vicinity of x — 0: 



u' (x = +0) - v! (x = -0) = -4v (x = 0) . 
v' (x = +0) - v (x = -0) = -4u (x = 0) 

Exact soliton solutions to Eqs. (ITU1) are looked for as 



{u(x), v(x)} = v^sech (s/2k {\x\ + {£, 77})) 



(11) 



(12) 



with £,?7 > 0. Symmetric and asymmetric solutions correspond to £ — rj and £ ^ ?y, respectively (the present model 
does not admit antisymmetric solutions). Total power ([7]) of solutions (|T2"T) is 

P = 2V2k 2 - tanh (V2kA - tanh (V2kr^j . (13) 

The substitution of expressions (|T2l into Eqs. (jlip yields the following equations which determine positive constants 
£ and rj: 



- sinh (V2k {£, r?}) sech 2 ( V2k {£, 77}) = sech (v 7 ^ {77, £} 



(14) 
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Using notation t\ t 2 = tanh (V2fcU,r?}), Eqs. (HI can be transformed into the following form: 

(ti - h) (t? + t\ + ht 2 - 1) = 0, (15) 
tita - 2/k, (16) 

where physical solutions are constrained to < t\ t 2 < 1. Equation (|15j) gives two solutions: ti = £2, which corresponds 
to symmetric solitons, and the other solution, which corresponds to asymmetric ones: 

tl + 4+ht 2 = l. (17) 

As it follows from Eq. ([To]) , the symmetric solutions have 

tanh (V2k {£, ?y}) = y/2/k, (18) 

hence they exist only for fc > 2, the total power (fl~3| of the symmetric soliton being 

^symrn = ^ ( Vk - Vt) . (19) 



The phase transition (SBB) occurs when the symmetric solution, t\ = ti = y/2~Jk, is simultaneously a solution to 
Eq. ||I7|), which yields 

k = 6, tanh (2y/3 {£, ??}) = -^=, P symm = (20) 
The solution for the asymmetric solitons, following from Eqs. (fT7)) and (fT6|) . is 



,,//-,, .\ fc - 2 ± Vfc 2 - 4fc - 12 
tanh 2 (a/21 {£, 77}] = ^ , (21) 

which exists exactly at k > 6, cf. Eq. (|20p . Further, the total power (| X3[) of the asymmetric soliton is 



P — 2 

1 asymm ^ 



2V2T- f y fc - 2 + Vfc 2 - 4fc - 12 + y fc - 2 - Vfc 2 - 4fc - 12 

and the relative asymmetry of the soliton is 



(22) 



P2+P2 



± Vfc-2 + \/fc 2 -4fc-12- v 7 ^ - 2 - Vfc 2 - 4fc - 12 
2\/2fc - f Vfc - 2 + /fc 2 - 4fc - 12 + Vfc - 2 - Vfc 2 - 4fc - 12] 



Figure [6] plots asymmetry O versus -Pasymm, as obtained from Eqs. (|23|) and (|22|) . The figure confirms that the SBB 
is indeed of the supercritical type when the linear coupling is spatially localized. 

The model based on Eqs. (j9|) is solvable too in the case of the self-dcfocusing nonlinearity, corresponding to the 
opposite signs in front of the cubic terms. However, in that case the model admits solely symmetric solutions, in the 

form of u(x) = v(x) = \/2fc/sinh (^/2k (\x\ + , with £ determined by equation tanh (^/2k^j = y/k/2, which has 



solutions for fc < 2, cf. Eqs. (JT2J) and dE 



V. CONCLUSIONS 



The objective of this work is to study the SBBs (symmetry-breaking bifurcations) of solitons in two-component 
systems which, unlike previously studied models, include the spatial modulation of the linear-coupling strength. To 
this end, two photonic models were considered, namely, the inverted virtual PhC (photonic crystal), and the fused 
dual-core spatial coupler. The former system is built as the periodic distribution of the density of dopant atoms, 
activated by the EIT, which induces the linear mixing between the two probe fields. The periodic density modulation 
makes all parameters of the medium periodic functions of the coordinate. Disregarding the XPM terms, we have 
found that the type of the SBB changes from sub- to supercritical with the increase of the total power of the probe 
beams. In the model of the fused dual-core coupler, the solutions for the two-component solitons were obtained in 
the exact form, the corresponding SBB being supercritical. 

The work can be naturally extended in other directions, including the interplay with the spatial symmetry breaking, 
and the consideration of higher-order solitons. A challenging possibility is to develop a two-dimensional generalization 
of the system. 
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FIG. 6: (Color online) Asymmetry O versus total power P aS ymm, in the solvable model based on Eqs. ([9]). 
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